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Abstract 

We present an analysis of four sum rules, each based on chiral sym- 
metry and containing the difference Py{s) — pa('S) of isovector vector 
and axialvector spectral functions. Experimental data from tau lepton 
decay and electron-positron scattering identify the spectral functions 
over a limited kinematic domain. We summarize the status of the ex- 
isting database. However, a successful determination of the sum rules 
requires additional content, in the form of theoretical input. We show 
how chiral symmetry and the operator product expansion can be used 
to constrain the spectral functions in the low energy and the high en- 
ergy limits and proceed to perform a phcnomcnological test of the sum 
rules. 



1 Motivation 



Despite a concerted effort by physicists extending over many years, an under- 
standing of QCD from first principles continues to be elusive. Fortunately, 
data continue to appear which provide a rather direct probe of the inner 
workings of the strong interactions. 

A case in point involves semileptonic tau lepton decay and hadron pro- 
duction in e^e~ scattering. For both of these, multihadron states such as 
27r, Svr, . . . are excited from the operation of quark vector and axialvector 
currents on the QCD vacuum. In this paper, we shall restrict our discussion 
to the isospin currents, fl 

Va^ = <t^rq and Ai:,=q^rim , (1) 

where a = 1,2,3 and q = [u d). Such current-induced processes provide 
information about the current bilinears, 

(0|r(Kr(x)F,'^(0))|0) and (o|r (A(^(xM|;(0)) |0) . (2) 

It has been long recognized that these quantities appear in certain chiral 
sum rules. Since these sum rules follow rather directly from QCD and 
chiral symmetry, a test of their correctness is, in effect, an experimental 
check on the validity of QCD itself. 

Unfortunately, there are several formidable obstacles to a successful im- 
plementation of this procedure. For one, the sum rules encompass an infinite 
range of energy, whereas existing data covers a very modest range, s < rn^ 
for tau decay and s < 5 GeV^ in e"'"e~ scattering. Moreover, as we shall 
see there are uncertainties in existing data which future experimental work 
must clear up. 

We feel that such difficulties can be overcome. In this work, we shall 
argue that a combination of chiral symmetry and QCD sum rule methods 
constrain the low-energy and high-energy limits of the dispersion integrals 
and that data can be used to fill in much of the rest. The content of the 
paper is organized as follows. We begin by introducing the spectral func- 
tions and their sum rules in Sect. 2, review the state of existing data in 
Sect. 3 and then describe various theoretical constraints in Sect. 4. Armed 
with experimental data and theoretical constraints, we present details of a 
phenomenological analysis of the chiral sum rules in Sect. 5. The final part 
in Sect. 6 summarizes our findings. 
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2 Spectral Representations 

In this section, we shall work exclusively in the chiral world of massless u, d 
quarks. Here, the spin-0 axial contribution is given by the pion pole and 
the two-current time-ordered products can be expressed in terms of spin-one 
spectral functions /?v,a('S))'0' 



{0\T{V,^{x)VaO))\0) 

(27r)4 pi-s + ie 



i5ab I ds py{s) i-sg^'^ - d^dn I — ^ — (3) 



and 



{0\T {A^,{x)Am)m = -idabF^d'^d'' 



, ^4 p-ip-x 



(27r)4 p2_^ie 

+i6at / ds PA{s){-sg^^ - d^dn / — . (4) 

For completeness, we also give the corresponding relations involving non 
time-ordered products. 



1 

2^ 
1 

2^ 



d^x e^'^-^{0\V,^{x)V,-mO) = iSab Pv{q^) {q^q" - q^g^ 
d^xe^''-^{0\A^{x)A'^mO) = iSab[pA{q^){q^q''-q'9^n (5) 



It follows from chiral symmetry and the high energy behavior of QCD 
that the vector and axialvector spectral functions contribute to certain sum 
rules, 

Hds Pv(g)-M(^) ^ -4Zio = (2.73 ± 0.12) X 10-2 (6) 
Jo s 

X) 

ds {py{s) - pa{s)) = Fl (7) 

= (8.54 ± 0.06) X 10-^ GeV^ , 

ds s {py{s) - pa{s)) = , (8) 

dss\n(j^^ (pyis) - pa{s)) = -^(ml^-mlo) (9) 

= -(6.19 ± 0.03) X 10-^ GeV^ . 
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In the remainder of the paper, we shall refer to the above sum rules re- 
spectively as WO, Wl, W2 and W3. In the first one, WO, the quantity 
Lio is related to the renormalized coefficient L^iq{p) of an 0{E'^) operator 
appearing in the effective chiral lagrangian of QC-D.@ Although the value 
of l[q^(/u), which is measured in the radiative pion decay vr — > euj, refers to 
a renormalization scale /x, the quantity Lio is itself independent of fi. 



Zio = 4o(a^) + 



144 



In 



-(6.84 ±0.3) X 10" 



(10) 



The next two relations Jl^l and W2, are respectively the first and second 
Weinberg sum rules. 'S'H The final sum rule, 1^3, is a formula for the vr^- 
7r° mass splitting in the chiral limit. M Although apparently containing an 
arbitrary energy scale A, this sum rule is actually independent of A by virtue 
of W2. For reference, we have displayed the physical values of the nonzero 
entries which appear on the right hand side of the chiral sum rules. These 
quantities have slightly shifted values in a chiral invariant world. This point 
is discussed at the end of Sect. 3. 

We note in passing that the current correlators defined in Eqs. (^Q) 
have been the subject of much recent attention. Several analyses have been 
carried out of hadron production in tau-lepton decay in order to obtain a de- 
termination of as(mT-), the running strong fine structure constant evaluated 
at the tau mass scale.! 



3 Data Inputs 

It is possible, in principle, to analyze the chiral sum rules on the basis 
of pure theory. For example, in the original derivation of W3, the pion 
electromagnetic mass difference was estimated by using p(770) and ai(llOO) 
contributions to saturate the vector and axialvector spectral functions. 0. 
However, a more sound procedure is to use data from tau lepton semileptonic 
decays into pions and/or pion production in e~^e~ annihilations. 

The rate for tau decay into an even or odd number of pions at invariant 
squared-energy s is given by @ 



ds Snml 



{m^ — s) 




ill) 
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and the corresponding nvr branching ratio is 



Bn-K = In-K With Imr= ds [ I 5- iH 5- p„7r(s) • 

Svrr^^all J{nm^f V J \ J 

In addition to the Particle Data Group (hereafter PDGllLiJ, the primary 
sources for tau decay data are the ARGUS and CLEG O detectors, 
via the reaction e^e~ — > T'''r~, and the LEP detectors il3 via the decay 
—>■ T~^T .A. review of tau physics up to 1988 is given by Barish and 
Stroynowski.lc3 Experimental aspects of tau decay continue to be presented 



up to the most recent conferences. 113' 113 112 

At first, using e~^e~ annihilation data to test the chiral sum rules would 
appear to be problematic. Although the range of energy is (at least in 
principle) unlimited, the production mechanism involves the electromagnetic 
current and so is generally a mixture of the needed isospin 1 component 
and an unwanted isospin component. However, for final states which 
consist of an even number of pions it is only the isovector electromagnetic 
current which contributes. This is a consequence of the G-parity relation 
GriTT = {—)"', together with the property that any nvr final state produced by 
the action of the electromagnetic current on the vacuum must have charge 
conjugation G = —1. Since G = G{—)^, it follows that 1 = 1 for n even. 
Of course, the isospin components which are measured in e~^e~ scattering 
and tau decay are distinct, having 1^ = and = 1 — i2 respectively. The 
corresponding 2tt and 47r states are related by 



r_|7r+7r") = V2\ 



o„-\ 



vr vr 



T_|2vr+2vr") = 2\/2|vr°7r+2vr") (13) 
r_|vr+vr"2vr°) = \/2|vr~3vr°) + 2\/2|vr+vr°2vr-) . 

Extraction of the nvr component of py{s) from e+e~ data proceeds via the 
relation 

For the specific case of two-pion production, the e~^e~ annihilation data 
is often expressed in terms of the pion electromagnetic form factor F'^'^{s) 
evaluated at squared-energy s. In this notation, one has 

Pr(^) = 4ilfl-^)'''|^'^^(^)l^ • (15) 
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A plot of the pion form factor in the timelike region appears in Fig. 1. Early 
results on pion production in e^e~ scattering is summarized in Ref. fl^ . 
Expemnente at Frascati, Orsay and Novosibirsk have continued to supply 
data.ii-il 

As a whole, tau and e^e~ data reveal that each of the multipion con- 
tributions rises fairly sharply from threshold to a peak value and then falls 
rather more slowly. At energies below 2 GeV, the role of meson resonances 
is significant. Thus, the 2ti contribution has the familiar narrow resonant 
structure of p(770), the Svr modes are dominated by ai(1260), and the 47r 
sector is influenced by /9(1450) and p(1700). Although lacking a detailed dy- 
namical understanding of higher multiplicity distributions, we can anticipate 
their form as a consequence of general physical considerations. They will 
occur sequentially in mass (due to increasingly higher mass thresholds), ex- 
hibit slowly decreasing peak values (from competition with other channels 
as constrained by unitarity), and become increasingly broad (since phase 
space grows with particle number). As energy increases, one soon enters the 
asymptotic domain and the sum over all modes becomes featureless, in a 
manner similar to the total e^e~ hadronic cross section. 

Let us now present a critique of the current status for various individual 
contributions: 

Tau Lepton Properties 

The most recent complilation given by PDG for the primary tau lepton 
properties of mass (mr), lifetime (r,-) and electron branching ratio = 

Tr^et'T-i'e /-^T— tall '^^^ 

^PDG _ ^ jg^ Q^y ^ ^PDG _ ^ ^PDG _ Qjygg ^jg) 

However, there have been recent downward revisions to'Ell 

rrir ~ 1.777 GeV , ~ 297 f s , ~ 0.1771 , (17) 

and these are the values that we shall use throughout our analysis. Observe 
that there is a slight inconsistency between the listed central values for m,-, 
Tr and Be in Eq. (p!2|). As a result, the theoretical constraint 

= "-"''''whir, (18) 



5 



is not exactly satisfied. Thus, the relation between nvr branching branching 
ratios and spectral functions in Eq. ( |T2|) depends on which of the following 
relations is assumed, 

(i) Bn^ = 4.39—4^^ In. or (ii) B^. = 71.31 B, . (19) 

lO^^'^ sec 

Throughout this paper, we shall for definiteness assume that i?e = 0.1771 
in our numerical work. This value is an average of the value in Eq. ( [T^ ) and 



the recent experimental determination cited in Ref. [28|. 

Two-pion Component 

Data for py" comes from the tt'vt'^ part of one^prong r decayEll EH and 
from the vr+vr^ final state in e^e~ scattering. '1^' EH The consistency 

of tau decay and e^e~ annihilation results in the vicinity of the p(776) peak 
has been verified by Gan in Ref. [^ ]. We display in Fig. 2 the two-pion 
spectral function p'y{s) as inferred from the pion form factor data of Fig. 1. 
Numerical integration of yields a 2it branching ratio of B2-K = 0.247. 
There are also recent determinations of the h^i r^ ( where /i is a hadron) and 
vr-vrO branching ratios by CL^olH and L^pH, 

^ Q 2483 ± 0.0015 ± 0.0053 , 
B^ePq = 0.243 ± 0.008 , (20) 
^ 0.2435 ± 0.0055 , 

where the t[~-k^ value is inferred by subtracting off the K*'^ branching ratio 
from that of h^n^. The above are in reasonable accord with the value cited 
by the PDG& , which is based on earlier data. 

All in all, the two-pion part of the vector current spectral function is 
well determined. As we shall now see, although much is known about the 
three-pion and four-pion distributions, more experimental input would be 
welcome. 

Three-pion Component 

We denote branching ratios for the two 37r modes in r decay as B^+2.,r- 



and S^-27r0- Basic isospin considerations imply the inequalities e1 

1<5^^4 l^E^^l . (21) 

2-^3^-5 5- B3. - 2 ^ ' 
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PDG lists the three-hadron branching ratios based on a number of experi- 
ments as 



Bh+2h- = 0.084 ± 0.004 and 3^-2^^ = 0.103 ± 0.009 . (22) 
CLEO^^ and LEP^^ have recently announced the results 



^ Q_Qg2i -t 0.0015 ± 0.0038 ± 0.0028 , 
Bh%o = 0.104 ±0.008, (23) 
Bj^'^^f^^ = 0.0949 ± 0.0036 ± 0.0063 . 

Collectively, these imply that B^+2t^- — B^-2-,^0 and indicate a total 37r 
branching ratio of 16%-19%. For definiteness, we shall assume that both 
modes have equal decay rates and employ total three-pion branching ratios 
in this same range. This is in accord with the findings of Davier who, in a 
review of early and recent 37r branching ratio determinations, summarizes 
the current situation aslI3 

B^+2^- = B^-2^0 = 0.0903 ± 0.0036 , (24) 

but at the same time makes the cautionary remark that 'extreme care should 
be exercized when using world average values for branching ratios'. 

The equality of r^+27r- and is expected of a final state which is 

dominated by the Ai resonance. The isospin decomposition 

\A,) = ^\p\-)-^\p-n^) (25) 

shows that the 7r"'"27r~ (from the first term) and the 7r~27r^ (from the second 
term) final states will occur with equal probability. 

As regards the spectral function p^, reconstruction from experiment 
would require Stt mass distributions for both tt'^Itt" and tt~27:^ modes. 
None of the latter exists. However, in view of Ai dominance it suffices 
to know just the 7r"'"27r~ spectrum. We refer the reader to Ref. ||l^ for 
histograms of the 3tt mass distribution measured some time ago by the 
DELCO, MAC and MARK II detectors. The literature also contains an 



early ARGUS determinationllll , corresponding to the rather small branch- 
ing ratio i?^4- oj— = 0.056 it 0.007. In this paper, we shall employ recent 
ARGUS datalMy to determine . From the number of counts AA^ per 
energy bin AE, one can construct the 3tt spectral function via 

a^^(s) = B,^ 1 1_AN 

'^^ ^ ' 24tt^V^^ Be (1 - s/mlf (1 + 2s/m2) 2^iVtot AE ' ^ ' 



7 



where E = ^/s. Upon taking B^^^ = 0.17 rather than the value given in 
Ref. (which would imply B^^^ = 0.13), we obtain the spectral function 
shown in Fig. 3. The large error bars near the end-point occur because one 
must divide the number of counts per energy bin by a phase space factor 
which vanishes at s = m^. Even so, it is clear from Fig. 3 that tau decay 
data is able to cover essentially all the region where is nonvanishing. 

Four-pion Component 

In tau decay, there are two Air modes, 7r~37r'' and 7r+7r*'27r^. The corre- 
sponding four-pion final states in e'^e~ scattering are 27r^27r~ and 7r~''7r~27r''. 
The four-pion spectral function measured in tau decay can be decomposed 
as 

Pv{s)=pV''i^) + PY'''i^) ' (27) 
where the quantities on the right-hand side are inferred from the four-pion 
mass distributions in the 7r~37r'' and tt'^tt^2tt~ modes. It is also possible to 
obtain the quantities py ^ and Py^^^ from e~^e~ — > Att cross sections via 
the relations 

Py^'i') = :iT^^^2.+2.- , (28) 

Py~~°i^) = 32^3q,2 ■^(^27r+27r- + 2tT^+^-2,rO ) • (29) 

The set of 7r~37r'' and > Sir^ tau decay branching ratios taken from 
recent conference presentations Refs. [|l^], [13] and the average cited by the 
PDG& provide a reasonably consistent picture, 

^ Q QQgg ^ Q_Qoo7 -t 0.0012 ± 0.0003 , 

^LEPg^o = 0.0153 ± 0.004 ± 0.006 , (30) 

^h"->3^o = 0-027 ± 0.009 . 

implying 5^-3^0 ~ 0.01. To our knowledge, there is no published spectral 
information from tau decay data for the r — > ■k~ Jit^^ Vj- mode. However, we 
can obtain Py^^^is) from (y2-K+2-K- ^ ™ Eq. (|28|). The set of 27r+27r~ cross 



section data taken from Ref. (for ^/s < 1.4 GeV) and from Refs. [20|, 
(for -^s > 1.4 GeV) is displayed in Fig. 4. Note that Fig. 4 clearly demon- 
strates how py" must extend beyond s = m^, i.e. tau decay alone cannot 
determine the Att spectral function over its full range. The resulting py^^^ 
obtained in this manner is shown in Fig. 5. From this, we predict the 7r~37r'^ 
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mass spectrum to have the form shown in Fig. 6, where a smooth curve has 
been added to help guide the eye. Integrating this mass spectrum leads to 
a branching ratio B^-^^o — 0.009, in good agreement with the values of 
Eq. (IgD. 

Turning to the larger 7r+27r~7r'' tau decay mode, we cite the branching 



ratio values appearing in Refs. |13[, [|10|], |lf 

B\^lh->i.^ = 0.0489 ±0.008, 

B^^^^-lo = 0.054 ± 0.004 ± 0.005 , (31) 
^r+°2h->i.o = 0.053 ±0.004, 

where the ARGUS value is preliminary. Two of these values include un- 
known amounts of strange particle contributions. Allowing for such non-47r 
contributions, the PDG and LEP values indicate a 7r+27r~7r'^ branching ra- 
tio in the 0.040 — 0.050 range. We are aware of just one published tt^2tt^tt^ 
mass spectrum measurement, an ARGUS analysis (Ref. [^^) which cites the 
branching ratio i?^+2^-^o = 0.042±0.005±0.009. This value is smaller than, 



although not inconsistent with, the more recent ARGUS value of Eq. (31). 
Using an appropriate binning procedure, we can construct the full An spec- 
tral function Py^{s) over the restricted energy region s < by combining 
the 7r+27r~7r'' tau data together with the 27r+27r~ cross sections. The result 
is shown in Fig. 7, together with an asymmetric Breit-Wigner fitting curve 
(see Sect. 4). 

Alternatively, one could use the combination of cross sections as in 
Eq. (|29| ) to determine the total 47r spectral function. In principle at least, 
this procedure can provide the 47r spectral function over a larger energy in- 
terval. 7r^7r^27r^ cross section data taken from Ref. p5| (for ^/s < 1.4 GeV) 
and from Ref. [^] (for y/s > 1.4 GeV) are displayed in Fig. 8. These turn 
out to imply a 7r"'"27r~7r'^ tau decay branching ratio somewhat smaller than 



the recent determinations cited in Eq. 31. For this reason, we have chosen to 
base our analytical work on the determination of Py^{s) as shown in Fig. 7. 

Higher Components 

The five-pion component in tau decay involves the branching ratios 
B^-^^o, B-^^~2tt+-> ^iid 527r-7r+27rO- Isospiu Constraints for these modes are 

- ^5^ " 10 ' - B^^ - 35 ' 35 - B^^ - ■ 

(32) 
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At present, Refs. |12], [10|, provide the following branching ratio deter- 
minations, 



0.0015 ± 0.0004 ± 0.0005 ± 0.0001 , 
0.00056 ± 0.00016 , (33) 
0.0054 ± 0.0005 ± 0.0008 . 

Noting that B^^~2tt+ ^ -^3/1-2/1+5 find that these values are consistent 
with the bounds of Eq. (|3^). However, the isospin bounds do not imply any 
useful information for disentangling -637^-277+ from -63/8-2/1+- To our knowl- 
edge, the above branching ratios are the only Stt data currently available. 
To obtain useful spectral infomration for the 57r mode requires a substantial 
number of events, e.g. as would be generated from a r factory. 

Some Ovr spectral information is available from e^e~ scattering and some 
Gtt branching ratio information is available, but in view of the paucity of 5tt 
data we have not included this sector in the analysis described in this paper. 

Although our statement of the chiral sum rules in Sect. 2 refers to the 
limit of massless quarks, the data reviewed in this section are taken in the 
real world of mu,d 7^ 0. In principle, we might attempt to perform corrections 
on the data set with an eye towards working in the chiral limit. For example, 
it is apparent that taking m^r — > would induce minor shifts in resonance 
masses and phase space. However, we anticipate that such effects would be 
of order /A^ with A ~ 1 GeV. Since such changes are much smaller than 
the present uncertainty in the data, it seems most prudent not to try to 
model the effect of finite mu,d effects. Thus, we shall use the unmodified 
experimental information in our phenomenological analysis. 

4 Theoretical Constraints 

Although the vector and axialvector spectral functions are not theoretically 
prescribed for all values of s, it is possible to place constraints on their low 
and high energy limits. 

We begin by taking the Fourier transform of Eq. (^), 

= Sab [q^iqu - 9f,uq^) nv(g^) . (34) 
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rCLEO 
-°h-47r0 
nPDG 

rARGUS 
^2tt-tt+2tt0 



and of Eq. 



(fx exp*5-^(0|T(A^(x)A^(0))|0) 



= 5ab {q^Qu - 9^.q^) HaC?') - q^q'^uf^ (q^) . (35) 

In the chiral limit, the spin-0 axial contribution I[^^\q'^) is given entirely 
by the pion pole. The correlators Uy^Aiq'^) are the real parts of analytic 
functions whose imaginary parts are the spectral functions /?v,a('S). Of in- 
terest to us are the dispersion relations involving the difference of vector and 
axial vector quantities, 



+ 



I ds s 



s — — le 
Pvjs) - Pa{s) 
s — q"^ — ie 



(36) 



In our normalization, the behavior of the individual fly, a (9^) to leading 
order at large q^ is 




(37) 



where // is the renormalization scale. In order to determine the difference 
nv(9^) — nA(g'^) for q^ large but finite, one must go beyond the form in 
Eq. (^7|). From the operator product expansion of vector and axialvector 
currents, we learn that the asymptotic dependence is 0{q~^) and the local 
operators which control this behavior are the four-quark condensates. 0^''^ 
In the approximation of vacuum saturation, one finds 



nv(g' 



-nA(^^) - 



9 



1 + 



q" 



4tt 



247 



+ In 



P 



(38) 



It turns out that theory also predicts the low energy or threshold behav- 
ior of the vector and axialvector correlators. However, in this case it is the 
machinary of chiral perturbation theory that is invoked to showE^' 



q 



2 + 48^2 



nv(g2) - IlAiq^ 
where L^^q (/i) is defined in Eq. ( p^ 



In 



P 



+ 



4^{P) 



(39) 
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The above statements all involve the correlators nv,A- Similar threshold 
and asymptotic constraints can be placed on the spectral functions py^A{s)- 
Thus, as stated in Ref. |Q the threshold behavior (s 4m^) of py{s) and 
Pa(s) is 

Pvis) ~ -±^(^1-^) e{s-Ami) + Oip') , (40) 
PAis) ~ 0(/) . (41) 

Later in this section, we shall sharpen the threshold result for pA by speci- 
fying the Svr threshold contribution in more detail. The perturbative result 
for the asymptotic limit s ^ oo of the individual /9v,A to leading order is 

/'v,a(^)-^(1 + ^) • (42) 
Finally, Eq. ( |38| ) determines the asymptotic form of Py{s) — PAis) to be 

Py(s) - pa(s) ~ ^ ^ ^^^3 for large s . (43) 

Note that the difference of the spectral functions is of order a^. Our analysis 
of the chiral sum rules will require a numerical value for the coefficient C of 
the term. From the estimate 

{y/a;qq)l ^ (0-24 GeV)*^ ~ 1.9 x 10"^ GeV^ (44) 

and taking ~ 0.2, we obtain 

C ~ 3.4 X 10"^ GeV^ . (45) 

The magnitude of this quantity is obviously model dependent and quite pos- 
sibly will be modified by future work. However, even folding in uncertainties 
of the vacuum condensates, it is clear that the coefficient of the term is 
very small, and that pv-PA approaches zero very quickly at large s. Even if 
one imagines pv{s)-pA{s) to exhibit increasingly damped oscillations indefi- 
nitely, duality suggests that Eq. (^3|) captures the correct average behavior. 
As we shall see, the asymptotic constraint of Eq. ( |4^ ) will have significant 
impact in the analysis of the chiral sum rules to come. 

The results presented in this section are well known, and collectively 
they represent a fairly powerful set of conditions regarding how the chiral 
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correlators can behave. Actually, additional thought can reveal even more. 
For example, let us employ the asymptotic behavior of Eq. (^) in the dis- 
persion relation of Eq. (|3^). Expansion of the dispersion integral in powers 
of yields the sum rule 



327r 
~9~ 



ds (pv(s) - Pa(s)) 



(46) 



Some care must be taken to interpret this result correctly. Observe that 
Eq. (H) is valid to 0{as). Since the 0{s~^) tail of py{s) - Pa{s) is itself of 
higher order in a^, one must subtract it off in this relation. Accordingly, we 
write 

327]- /"oo 

— {^/o^sqq)l = - ds {py{s) - Pa{s))' , (47) 

y JO 

where {py{s) — PA^s))' refers to subtraction procedure just mentioned. 

Finally, let us return now to the matter of the the threshold behavior 
of /5a(s)- It is associated with the Svr component. Using chiral lagrangian 
methodsfl, we have determined that in the chiral limit the threshold be- 
havior is 



(0|A'^(0)K+^vrO„V 



(p^+P^-2p^ 



3sj 



(48) 



where Q = pj^ +P- +Po-, s = and 
of terms, this can be expressed as 



3F^s 

{p+ + P- ) ^ • Up on rearrangement 



(0|A'^(0)|7r, 



2 



Q-Po 



Po 



(49) 



Observe that this obeys {0\d^A^^{0)\37r) = 0, as must be the case since it 
is only the spin-one part of the axial current which contributes here. The 
threshold behavior of p^{s) can then be read off from the general form of 
Eq. (^) by first squaring Eq. ( ^9|) and integrating over 3tt phase space. 



Pfis) 



967r(47rF^)2 



+ 



(50) 



In practice, however, the very low s behavior for the Svr state turns out to 
be less important in the chiral sum rules than the higher energy effect of the 
Al resonance. 
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5 Empirical Determinations of Chiral Sum Rules 



Before discussing our own methodology, we wish to take note of two inter- 
esting works involving aspects of chiral sum rules. In the earlier of these'0, 
spectral functions based on 27r, Svr and 47r data extracted from tau decay 
are used to study the sum rules Wl, W2 and W3. The associated spectral 
integrals are studied as a function of cutoff sq for sq < 2.5 GeV^. However, 
the Svr analysis was taken from Ref. |11] data, which as we have seen is not 



in agreement with a large number of more recent branching ratio measure- 
ments. In addition, the 27r spectral function was too small in the vicinity 
of the /o(770) peak by roughly a factor of 2. Interestingly, these two fea- 
tures combined to make the sum rules appear reasonably in agreement with 
expectations, although this result was fortuitous. 

In the analysis of Ref. [^] , 2-it and Air data are inferred from e'^e~ scatter- 
ing, but the Stt data again comes from Ref. In fitting these components, 
use is made of Breit-Wigner resonance forms with energy-dependent decay 
widths. The particular form of energy dependence is taken from Lorentz- 
invariant phase space. The Svr and higher components are parameterized to 
give rise to the asymptotic behavior 

5 

Pv,a(s) ~ and pv - M ~ ^(s^^) • (51) 

This disagrees with the chiral and operator product expansion results of 
Eq. (42) and Eq. ( |43| ) respectively. The specific form in Ref. [|l^] used for 



the higher components is chosen to fit the sum rule Wl exactly and the sum 
rule WO is evaluated in terms of the fit. 

The summary given in Sect. Ill of available data demonstrates that there 
is hope for successfully extracting much about the spectral functions py^a{s) 
from experiment, but that our knowledge of them will always be limited. In 
view of the current database, we have decided that for the purpose of testing 
the chiral sum rules, it is most prudent to take the empirical 27r, Svr and 47r 
modes explicitly into account, and to treat all higher components according 
to some reasonable prescription. We have followed two distinct approaches 
in doing so: 

1. Numerical: We ensure that the empirical databases for pv,a('S) evolve 
smoothly in the variable s to the correct asymptotic limits by gener- 
ating smooth curves separately for py and pA which pass through the 
experimental data sets at low energy and which satisfy the asymptotic 
limits, while reproducing the four chiral sum rules. 
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2. Analytical: This actually encompasses a class of fits to the differ- 
ence of spectral functions py — p\ in which all contributions higher 
than the four-pion sector are lumped into a single theoretical term. 
A convenient method for constructing the spectral functions this way 
is to start with a delta function form, then introduce finite-widths 
via Breit-Wigner representations, and finally modify these to a more 
realistic asymmetric form. 

Let us consider each possibility in turn. 

Numerical Representation 

Surely the simplest method for generating acceptable global versions of 
pv,A is to numerically smoothly join the low energy empirical data with 
asymptotic theoretical information. A reasonable region for matching the 
two occurs at about s ~ 4 — 5 GeV^. As we have already seen, the constraint 
of Eq. ( ]4^ ) reveals that the spectral functions approach a nonzero constant 
at infinite energy, with an additive correction factor proportional to as{s). 
In general, we expect the large s behavior 

V,A V,A 

pyAs)=p'o'^ + ^ + ^ + ... , (52) 

We can obtain a determination of the dominant power correction in the large 
s limit as follows. Using an operator product expansion, Braaten, Narison 
and Pich have displayed the structure of the correlators nv,A for Euclidean 
momenta —Q^ in Appendix A of their paperlH, 

n r n2^ - ^v,a.^x Oi'^iQ) (h'^iQ) 

iiv,A(-V ) - % [Q) H Q2 ^ gi ^ • • • > (53) 

where s = —Cf' and each of the a^'^{Q) is expandable in powers of as{Q)- 
We require the imaginary part of this expression, analytically continued to 
timelike momenta. Since the 0{Q~'^) contribution to nv,A turns out to be 
proportional to quark mass, the coefficient p^'^{s) vanishes in the chiral 
limit. Thus, the first nonvanishing subleading contribution is the 0{Q~^) 
component, from which we extract the result 

p2'^is) ~ ^^^^^ ^ ■ 



15 



Let us estimate the magnitude of this quantity at s = 5 GeV^. The 
strong fine structure constant is determined in terms of an assumed value 
for the QCD scale parameter, 

Aqcd = 150 MeV a^i^ GeV^) ~ 0.26 (55) 

and the gluon condensate from phenomenological applications of QCD sum 
rules, 

(^G2^o = (0.02 ± 0.01) GeV^ . (56) 

TT 

Altogether, these values imply that the 0{s~'^) component to pv,A has a 
tiny coefficient, 

p^'^ ~ 0.9 X 10-^ GeV^ . (57) 

As a result, even at the modest energy s ~ 5 GeV^ the 0{s~'^) term has a 
negligible effect. Of course, this common addition to pv and p\ will cancel 
when the difference is taken and hence will not contribute to the chiral sum 
rules. In addition, the difference of the spectral functions was chosen to 
be compatible with the asymptotic constraint given in the previous section. 
Again the magnitude of this contribution is so small that it is essentially 
irrelevant at modest energies. Thus while we have made an effort to generate 
spectral functions with the right high energy behavior, the precise value of 
the high energy terms is not important since their numerical size is quite 
small. 

Let us summarize our procedure at this point. We have generated nu- 
merical representations for pv and pA which fit all available data on mul- 
tipion production, and which are compatible with the theoretical behavior 
expected at high energy, and which when integrated yield the correct exper- 
imental values for the four chiral sum rules. Although highly constrained at 
low and high energies, the spectral functions have a modest uncertainty in 
the s = 2 — 4GeV^ range, and this was exploited in order to precisely dupli- 
cate the expected values of the sum rules WO, Wl, W2, W3. Of course, both 
this numerical procedure and the analytic one to follow arc subject to the 
choice of input values. Given the uncertainties in especially the Stt and 47r 
branching ratios, we have explored a variety of possibilities. Corresponding 
to the input set 

= 0.240 , Bs^ = 0.165 , B^^ = 0.048 , (58) 
curves for pv, PA and pv-PA are displayed respectively in Figs. 9-11. 
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Analytical Representations 



(a) Delta Function 

In the delta function description, the spectral functions are 

py{s) - pa{s) = X: {-)'^'Fi6{s - ml) . (59) 

k=l 

This representation of the spectral functions, although crude, provides a 
nice pedagogical example with which to organize one's thoughts. We begin 
by noting that there are 2 parameters per contribution, a mass m and a 
coupling F. In this approach, the four sum rules reduce to 

WO = ' (60) 



k=l 
4 



mt 



Wl = i-f-'Fl = Fl , (61) 

fc=i 

W2 = i-f^'^l Fi = 0, (62) 



k=l 
4 



W3 = Y i-f^'^l Fk Inmi , (63) 



k=l 

and the tau decay branching ratios become 

^ rrifX ^ ( ^ 2ml 



Bk = 71.62 B,F^ [1-^] 1 + ^ (A: = l,2,3). (64) 

m£ / \ mt I 



We can specify the mass parameters from the observed nvr {n = 2, 3, 4) 

distributions and the coupling strengths -^1,2,3 from the observed branching 
ratios (c/. Eq. ([l2|)). Even in the extreme narrow width approximation of 
the delta function representation, this step turns out to be a surprisingly 
accurate one, e.g. finite width effects modify the branching ratio relations 
by only a few percent. 

This leaves the problem of determining the coupling F4 and mass m4. 
It seems reasonable to require that the sum rules Wl and W2 be obeyed 
exactly (c/. Eq. (|6l|) and Eq. (|62[)). Thus we determine F4 from Eq. (|6l|) 
and m4 from Eq. (|6^). Doing so leaves the remaining sum rules WO and 
VF3 as predictions. 
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(b) Breit-Wigner 

In the finite width Breit-Wigner (BW hereafter) extension of Eq. (|5' 
the spectral functions are represented by 



pvis)- pAis) 



4 

E 

k=l 



\k+l 



vr 



ml 



)2 + {mkVk) 



(65) 



Each contribution still has mass and coupling parameters and F^, but 
now in addition a width Ffc. In the most general Breit-Wigner resonance 
form, decay widths are taken to be energy-dependent, e.g. the 27r widths 
have the 0{p^) threshold behavior to reflect the P-wave nature of the 27r 
system. However, in order to maintain a simplicity of description and the 
ability to represent each rnr component in analytic form, we shall employ 
energy-independent widths throughout. 

An advantage of the delta-function approach was the ability to express 
all the sum rules and the tau branching ratios in elementary form. This 
is still true in the BW approximation for the sum rules WO, Wl, W2 and 
the tau branching ratios. For example, working with a generic Breit-Wigner 
spectral function. 



Pbsn{s) 



(66) 



where r = r/(7r?n,) is an expansion parameter for finite-width effects, we can 
evaluate integrals such as 



r(0) 
'BW 



ds /9Bw(s) 



F 



1 vrr 

1 arctan — -- — ^ 

TT {Km?) 



1 vrr 

— arctan — - — ^r— 

TT 1 — (So/rra^j 



,(67) 



and 



f(i) 
'bw 



ds spBwis) 



m^4l + F\ 



A /(l-mVA)VmVrVA2^ 
In — + 0.51n ^ ' ' \ 

m \ (1 - so/m^)^ + vr^r^ 



(68) 



By passing to the limit r = of zero decay width, we regain the results of 
the delta function model. 
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Analogous to the procedure used in the delta function model, parameters 
for the first 3 BW poles are obtained by fitting to experimental data. In 
passing, we note that since a given spectral contribution has the asymptotic 
behavior 0{s~'^), the individual BW spectral integrals for W2 and W3 are 
divergent. The numerator of the fourth pole is fixed by demanding that the 
C(s~^) asymptotic term in py-pA vanish. This implies 

^ {-)'+'FimkT, = . (69) 
fc=i 

By virtue of this relation, the asymptotic behavior of the vector and axi- 
alvector spectral functions becomes 

py(s)-pAis) = ^ + 0{s^) (70) 

where 

Cbw = -E H'^'FimlTk . (71) 

The constant Cbw ™ay be considered either as a quantity to be fixed by 
Eq. ( ^3|) or as a prediction of the analysis. 

We find it hard to see how any analytical study of the sum rules could 
succeed without incorporating the above features or something similar. In 
our approach, all the chiral sum rules are convergent even though individual 
pole contributions may diverge. Moreover, the condition given by Eq. (]69| ) 
ensures that our description has the smoothness in energy expected from 
duality. 

(c) Asymmetric Breit-Wigner 

Although having the virtue of simplicity, the representation of Eq. (|65| ) 
is deficient in several important respects. A Breit-Wigner form is symmet- 
ric about its resonant energy whereas the rnr contributions to pv,A exhibit 
asymmetric bumps. Besides, the BW contributions extend to energies lying 
below thresholds which characterize the various rnr components. 

An improved treatment can be realized in a variety of ways. One simple 
parameterization which treats the various components uniformly, yields a 
reasonable fit to data and allows us to maintain analytic control is 

^ ^ \^ r \fc+i Pk{s)F^mkTk r^^-. 
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The functions Pk{s) are polynomials 




2 \ 



(73) 



each containing two parameters, a threshold energy and an integer-valued 
exponent Uk = 1, 2, . . .. Although the choice = 1 is the simplest one, it 
yields an rnr spectral function which rises linearly just above threshold. This 
does not appear to provide an adequate fit to the data, and thus we have 
used nfc = 2 (fc = 1 . . . 4). A fit of this type to the Svr spectral function 
appears in Fig. 12. 

Introduction of the polynomials Pk{s) will clearly lead to integrals in- 
volving inverse moments of Breit-Wigner forms. 



Although it is straightforward to evaluate such integrals analytically, the re- 
sulting expressions can be quite cumbersome. In practice, it is more efficient 
to employ the recursion relation 



Upon applying relations of this type, we have carried out the calculational 
program described earlier in this section. As stated earlier, because of the 
uncertainty in the experimental Svr and 47r contributions, we have performed 
the analysis for several different sets of tau lepton branching ratios. Typical 
results are shown in Table 1, where we display both input values (parameters 
for the 27r, Svr and 47r asymmetric BW poles and the associated branching 
ratios) and results (parameters for the 4*'' pole, the coefficient C and the 
values for WO and VF3).'^ We have purposely exhibited two solutions (^^1 
and 7^2) to show that it is possible to fit the sum rules yet not obtain an 
acceptable value of C. The graph of py-pA corresponding to solution #3 of 
Table 1 is displayed in Fig. 13. The overall appearance of this curve clearly 
mimics the one in Fig. 11 which is based on the numerical approach. We 
view the fact that the different methods generate rather similar spectral 
functions as an indication that there is not great freedom in pv ~ PA once 
the theoretical and experimental constraints are imposed. 




(74) 




(75) 
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6 Concluding Remarks 



What we have done in this paper is to suggest a procedure for compar- 
ing four well-known sum rules of chiral symmetry with data from the real 
world of experiment. In addition, we have phenomenologically constructed 
the vector and axialvector spectral functions to the extent possible by using 
the full collection of available tau decay and e^e~ cross section data. To 
be specific, the data which constrains the spectral functions includes the 
e+e~ — > TT+vr", 27r+27r~ and 7r+7r~27r'' cross sections, the r 27r, Stt and 
47r branching ratios, and the energy spectra for Tr^Tr*^, 27r~7r+ and 27r~7r+7r'^ 
final states in r decay. Theoretical constraints include chiral symmetry at 
low energy, isospin relations for handling the data and the operator product 
expansion of QCD at high energy. It is clear to us that this activity will be 
repeated by ourselves or by others in the future. That is, we (perhaps opti- 
mistically) anticipate the emergence of improved data which will provide a 
yet more reliable foundation upon which to base the phenomenology. How- 
ever, we feel that our results are the the best determination of the spectral 
functions that can be made at this time. 

As regards the data, we urge that efforts be made to improve the deter- 
minations of the 3tt and Att contributions to the spectral functions. The Sir 
component can only be inferred from the tau decay hadronic distribution. In 
this paper, we were fortunate to have access to the recent ARGUS determi- 
nation of . Cross checks are always welcome, so we urge that data from 
tau production at both LEP and CESR be analyzed to extract the spec- 
tral function p^. In addition, the need for an improved 47r determination 
in e+e~ cross section data is especially acute for the 7r'^7r~27r° final state. 
Additional information on the 47r tau decay modes (7r'^7r°27r~ or tt^Stt'^) 
would also be welcome. 

In the latter part of our paper, we addressed the question of whether 
experimental data is consistent with the chiral sum rules. On the basis 
of our study, we conclude that existing data is indeed consistent with the 
chiral sum rules. It is important to not misinterpret this remark. Of course, 
since physical data will always be less than perfect, it is not possible to claim 
'proof of validity for the set of sum rules. It is evident to us that to insist on 
such proof would be foolhardy. However, given the number of constraints 
on the spectral functions, it is far from trivial that all of the chiral sum 
rules can be satisfied. Agreement at the level we obtained is an affirmation 
of the subtle and complex theoretical intuition (involving chiral symmetry, 
dispersion relations, the operator product expansion and the asymptotic 
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behavior of QCD) that leads to the sum rules. 
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Trial 


#1 


#2 


#3 


Two-pion Inputs 


Ml 


0.763 


0.763 


0.763 


/i 


0.157 


0.157 


0.157 


Ti 


0.123 


0.123 


0.123 


Sl 


4ml 


4ml 


4ml 


B2tt 


0.243 


0.243 


0.243 


Three-pion Inputs 


M2 


1.117 


1.117 


1.117 




0.244 


0.234 


0.250 




0.470 


0.470 


0.470 




0.500 


0.510 


0.550 




0.185 


0.168 


0.176 


Four-pion Inputs 


M3 


1.500 


1.500 


1.490 


/3 


0.192 


0.188 


0.215 




0.564 


0.700 


0.765 


S3 


0.700 


0.710 


0.690 


Batt 


0.048 


0.041 


0.053 


Output Values 


M4 


2.288 


2.055 


1.869 




0.068 


0.083 


0.116 


r4 


0.221 


0.750 


0.873 


Bn>4n 


0.0001 


0.0015 


0.0046 


WO 


0.0261 


0.0266 


0.0268 


W3 


-0.0062 


-0.0062 


-0.0062 


C 


0.013 


0.003 


3.8 X 10-5 



Table I: Asymmetric Breit-Wigner Representation. 
The energy unit is GeV. 
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Figure Captions 



Fig 


. 1 


Timelike pion form factor 


Fig 


. 2 


27r vector spectral function 


Fig 


. 3 


Svr axialvector spectral function 


Fig 


. 4 


Cross section for e+e~ — > 27r+27r~ 


r ig 


. 


p-w inierreQ iroin e e scatTjeimg 


Fig 


. 6 


Predicted Att mass spectrum in r — 7r~37r°^'T- 


Fig 


. 7 


/jy" from tau decay and e+e^ scattering 


Fig 


. 8 


Cross section for e+e~ —>■ 7r"'"7r~27r° 


Fig 


. 9 


Numerical fit to py 


Fig 


. 10 


Numerical fit to p\ 


Fig 


. 11 


Numerical fit to py — pA 


Fig 


. 12 


Fit of asymmetric Breit-Wigner solution to Sir spectral function 


Fig 


. 13 


Fit of asymmetric Breit-Wigner solution #3 to py — pA 
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